Role of Antisymmetric Exchange in Selecting Magnetic Chirality in Ba3NbFe3Si20i4 
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We present electron spin resonance (ESR) investigation of the acentric Ba3NbFe3Si2 0i4, fea- 
turing a unique single-domain double-chiral magnetic ground state. Combining simulations of the 
ESR line-width anisotropy and the antiferromagnetic-resonance modes allows us to single out the 
Dzyaloshinsky-Moriya (DM) interaction as the leading magnetic anisotropy term. We demonstrate 
that the rather minute out-of-plane DM component dc = 45 mK is responsible for selecting a unique 
ground state, which endures thermal fluctuations up to astonishingly high temperatures. 

PACS numbers: 75.30.Gw, 76.30.-v, 76.50.+g 



Spin-orbit coupling generally has no relevant effect in 
ordinary magnets, except to pin the orientation of spins 
with respect to the background lattice in magnetically 
ordered phases, singled out by the isotropic Heisenberg 
exchange interactions. An exception is provided on ge- 
ometrically frustrated lattices where the isotropic inter- 
actions alone are often unable to raise the macroscopic 
degeneracy of a ground state (GS) leading to uncon- 
ventional cooperative electronic states [iHll- Then, the 
magnetic behavior might be exclusively driven by minute 
perturbing terms in the form of anisotropic interactions 
emanating from the spin-orbit coupling. For instance, the 
antisymmetric Dzyaloshinsky-Moriya (DM) exchange in- 
teraction Q can induce order on a kagome lattice in the 
classical limit and can lead to quantum criticality 

in the quantum limit Furthermore, DM is known to 
alter phase diagrams of frustrated ladders [lOi and trian- 
gular lattices [ll| . It is also responsible for spin chirality 
|12| that can be long-ranged even in the absence of a 
classical magnetic order [13| . Spin chirality is one of the 
key concepts in the physics of strongly correlated elec- 
trons, as it is related to various intriguing phenomena, 
like the realization of s pin liquids [14|, magnetic-order 
induced ferroelectricity [15|, anomalous Hall effe ct JlGll 
and possibly high-temperature superconductivity |17|. 

In this context, Ba3NbFe3Si2 0i4 (BNFSO) is ex- 
tremely appealing due to its remarkable magnetic prop- 
erties. It crystallizes in a non-centrosymmetric trigonal 
unit cell (P321 symmetry). The Fe^+ {S = 5/2) spins 
reside on vertices of equilateral triangles arranged into a 
two dimensional (2D) triangular lattice (crystallographic 
ab planes in Fig. [1]) [18j]. The dominant exchange in- 
teractions are antiferromagnetic - the Curie- Weiss tem- 
perature is 9 ^ —180 K [18|-l20j - and thus frustrated. 
Nevertheless, a long-range-ordered (LRO) state is real- 
ized below T/v = 26 K, characterized by a 120° spin ar- 
rangement on each triangle. The moments are bound 
to the ab planes and form a magnetic helix along the 



crystallographic c axis, corresponding to the magnetic 
propagation vector q = (0, 0, r), r '-^ 1/7 [l8| . 

The system has been drawing considerable attention 
because of its magnetoelectric and multiferroic proper- 
ties [19rL21i] . However, its uniqueness stems from its dis- 
tinctive magnetic order and magnetic excitations. The 
LRO magnetic GS is doubly chiral [Fig.[ljb)] and single- 
domain - a single triangular vector chirality eA together 
with a single helicity en is chosen from four possible 
states in a structurally chiral crystal - which is be- 
lieved to be crucial for its multiferroic properties 
Moreover, below T/y one of the two magnetic excitation 
branches emerging from the magnetic satellites is com- 
pletely chiral over the whole energy spectrum, implying 
an unprecedented absence of chirality mixing of the spin 
dynamics 
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The chiral correlations remain present far 



23j . The helicity of the GS can be rationalized 



within the isotropic Heisenberg model if three competing 



(a)3' 




^^^^ 



,1 c 



FIG. 1: (color online). 2D "triangular" arrangement of Fe'^"'" 
in BNFSO. (a), (c) Exchange interactions Ji_5 and magnetic 
anisotropics of the Dzyaloshinsky-Moriya (d) and single-ion 
(D, F) type. The latter vectors are parallel to a local two- 
fold rotational axis defining the Z axis, (b) Two possible 
double chiral ground states, with eA and en denoting trian- 
gular chirality and helicity (with ^ 27r/7 pitch along c axis), 
respectively. 
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iiiterplane interactions [J3-5 in Fig. [IJc)] are considered. 
However, the isotropic exchange stih allows 18| for two 
(e//7 ea) pairs [Fig.[TIb)], and fails to justify the selection 
of the ab easy plane. Thus, the fundamental question 
about the mechanism, responsible for the experimentally 
observed unique single-domain chiral GS [18| in BNFSO, 
arises. 

It was suggested that the DM anisotropy might be re- 
sponsible for selecting the chiral GS [18| and for opening 
a small gap in magnetic excitations [22] of BNFSO. In 
another inelastic neutron scattering (INS) study the lat- 
ter was proposed to arise from single-ion (SI) anisotropy 
[2^ . A clear and quantitative determination of the source 
of the magnetic anisotropy and of its impact on the GS 
is therefore needed. In this Letter, we present our elec- 
tron spin resonance (ESR) study, which provides a di- 
rect insight to magnetic anisotropy of BNFSO. Jointly 
simulating the experimental ESR line-width anisotropy 
in the paramagnetic phase and magnetic excitations in 
the ordered phase, we determine the dominant magnetic 
anisotropy, which is of the DM type and has sizable both 
the out-of-plane (dc) and the in-plane (dab) component 
(Fig. [1]). Our mean-field based calculations single out dc 
as being responsible for selecting the unique GS. 

ESR is an extremely powerful technique for quantify- 
ing magnetic anisotropy 01 • The isotropic Heisenberg 
exchange commutes with the Sz spin operator (z de- 
notes the quantization axis set by the applied magnetic 
field) and therefore leads to a i5-function resonance. Fi- 
nite magnetic anisotropy then yields a finite ESR line 
width. In single crystals the line-width anisotropy in 
the paramagnetic phase can often unveil the dominant 
magnetic anisotropy, since different types of anisotropy 
reflect local symmetries distinctively, as they arise from 
different microscopic origins. Finally, ESR can also de- 
tect collective magnon modes in a LRO state. We there- 
fore performed an extensive ESR investigation on high- 
quality BNFSO single crystals in the temperature range 
between 4 K and 500 K. X-band (at 9.3 GHz) spectra 
were recorded on a home-build double-cavity spectrome- 
ter, equipped with a helium- flow cooling and a preheated- 
nitrogen-fiow heating system. Measurements at frequen- 
cies between 50 GHz and 400 GHz were performed on a 
couple of custom-made transmission-type spectrometers 
at NHMFL, Florida. 

In X-band, Lorentzian-shaped spectra are observed 
down to T/v (Fig. [2]), in accord with strong exchange 
narrowing regularly encountered in dense magnetic sys- 
tems. Below T/v, the paramagnetic signal rapidly disap- 
pears. Our calibration of the ESR intensity, which is pro- 
portional to spin-only magnetic susceptibility xesr, and 
the scaling between xesr and the bulk susceptibility Xb 
[Fig.EIb)] prove that the ESR signal of BNFSO is intrin- 
sic. The observed small deviations of the g-factor from 
the free-electron value go = 2.0023 are typical for Fe'^+ 
with a nearly pure ^S5/2 orbital sing let GS 21. A uniaxial 
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FIG. 2: (color online). X-band ESR results: (a) Comparison 
of ESR, Xesr, and bulk, xt, susceptibility for B||c. Inset: 
room-temperature ESR spectrum (blue) and the correspond- 
ing Lorentzian fit (red), (b) Temperature dependence of line 
widths and the line width ratio (inset). Angular dependence 
of (c) the line width and (d) the g-factor measured at 500 K, 
with corresponding fits (solid lines) explained in the text. 



anisotropy of the g-tensor and the line width [Fig. ^c- 
d)\ are justified by a three-fold rotational symmetry of 
iron triangles. Our fit of the factor anisotropy (25| 
yields the three eigenvalues of the g-tensor: gxx = 2.000, 
gvY = 2.001 and gzz — 2.013, and thus reveals that also 
locally the symmetry of the g-tensor is very close to being 
uniaxial. The polar axis Z on each site is set by a local 
two-fold rotational axis lying within the crystallographic 
ab plane (a direction for site 2 in Fig.[T]). The other two 
principal axes X and Y for site 2 lie 30° and 120° from 
the b* axis in the crystallographic b*c plane [25)]. 

The anisotropy of the g-factor reflects mixing of ex- 
cited orbital states into the ground orbital singlet of 
the Fe'^+ ion, which is induced by a spin-orbit coupling. 
The anisotropy of the line width arises from magnetic 
anisotropy present in a spin Hamiltonian, which origi- 
nates from the same perturbation. We address the issue 
of the magnetic anisotropy in BNFSO in the framework 
of the spin Hamiltonian 
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s,^^lBB■g■J2^,+n', (1) 



where, the first sum runs oven the spin pairs connected by 
one of the five different exchange interactions (Fig.[T]) and 
represents the Heisenberg term He, the second sum gives 
the Zeeman coupling Tiz and the third sum the magnetic 
anisotropy 7i' . The two dominant contributions 26j to 
the latter for Fe'^^ are the single- ion anisotropy [2 
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FIG. 3: (color online), (a) Selected ESR spectra at 97.6 GHz. 
Inset shows two distinct resonance modes at 4 K. (b) Tem- 
perature dependence of the central field of the two observed 
modes. 



and the Dzyaloshinsky-Moriya anisotropy Q 
"H'dm = dij ■ Si X Sj. 



(3) 



The symmetry of the SI anisotropy is the same as 
of the (7-tensor [S^l, therefore we omit the additional 
E(Sj^-S^) term. Since Ji is dominant [H, [111, we 
consider the DM interaction only between nearest neigh- 
bors. The DM vector may possess a non-zero out-of- 
plane dc and in-plane dob component parallel to the bond 
[Fig[TJa)]. The third component is forbidden by the two- 
fold rotational axis passing through the middle of each 
bond. 

Using the Kubo-Tomita approach [1], we can analyti- 
cally calculate the line- width anisotropy ABoo{0) at in- 
finite temperature, separately for both types of mag- 
netic anisotropy [25|. At finite temperatures AB{d) — 
c{0,T)ABoo{O), where c{e,T) differs from unity due to 
finite spin correlations. In BNFSO the uncorrelated para- 
magnetic state is not reached yet even at 500 K (~ 39) 
[Fig. mb)], which complements specific heat measure- 
ments showing that 20% of magnetic entropy is still miss- 
ing at 200 K [20|. This needs to be contrasted with sev- 
eral Cu-based 2D frustrated lattices, where the ESR line 
width was found constant ioi T >6 2^, 30| . 

We find the fits of the experimental line width 
[Fig. Wis)] virtually indistinguishable for the SI and the 
DM model. At 500 K we extract ^/c\D\ = 0.53 K, a/D = 
1.06, F/D = 0.09 for the SI model and ^/c\dc\ = 0.13 K, 
dab/dc = 2.6 for the DM model. It is worth noting that 
since Fe^+ is in the orbital S state, in the SI model only 
the a term is allowed in purely cubic, octahedral or tetra- 
hedral environments. In BNFSO the local symmetry is 
lower, therefore, D and a can have similar magnitudes 
Although the absolute size or the sign of the anisotropy 
terms cannot be revealed from these fits due to unknown 



c > 1, the constant ratio Ai?_L/Ai3|| = 1.27, observed 
in BNFSO above ~ [inset in Fig. [2fb)], unambiguously 
sets the ratio of the anisotropy terms for both models 
[2^. The temperature-independent line- width ratio re- 
veals that correlations are isotropic in spin space. 

In our attempt to determine and quantify the dom- 
inant magnetic anisotropy in BNFSO, we now turn to 
high-frequency measurements. In contrast to X-band, 
we can follow the ESR signal across Tn [Fig. [3l^a)[. The 
resonance broadens and shifts to lower field by ^ 2 T 
below Tjv. The temperature dependence of the line shift 
mimics an order parameter [Fig. [3I^b)[, and corresponds 
to the opening of a zero-field gap in magnetic excitations 
(Fig. |4|). This transformation of the paramagnetic sig- 
nal into a collective antiferromagnetic resonance (AFMR) 
mode Si is accompanied by emergence of another AFMR 
mode (S2) in our frequency window [inset in Fig. [3l^a)[. 
The frequency- field diagram of the two modes is shown 
in Fig. m for both relevant directions of the applied field 
with respect to the crystallographic c axis. 

We have simulated the AFMR modes experimentally 
detected at 4 K in a mean- field approximation (25| , again 
for both types of magnetic anisotropy [Fig.[4l^a), (b)]. We 
note that the zero-field gap to the lowest gapped branch 

51 is opened by finite anisotropy {dc or D, a, F) in a 
square-root fashion [25| , while the gap to the next branch 

52 is present already for isotropic He and only slightly 
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FIG. 4: (color online). Frequency-field diagram of measured 
(symbols) and simulated (lines) AFMR modes for (a) DM and 
(b) SI anisotropy at 4 K. (c) Field-induced phase transition 
predicted by the SI anisotropy model at 5.5 T for B_Lc. (d) 
Angular dependence of the two resonance modes at 4 K; mea- 
sured (symbols) and simulated for the DM (solid lines) and 
SI anisotropy (dashed lines). 
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shifts with anisotropy. The latter gap is reproduced with 
the exchange parameters Ji = 11.9 K, J2 = 3.4K, J3 — 
0.72 K, J4 = 0.24 K, and J5 = 3.4 K, which corresponds 
to a 20% increase [31| of the parameters obtained in the 
INS study 

The SI model [Eq. ([S])], at first sight, yields a rea- 
sonable description of the experimental v — B diagram 
[Fig. Hlb)] for D ^ -17 mK, a = -18 mK and F = 
—0.9 mK ■ However, a closer look at this model un- 
veils its inconsistencies with the experiment. The simu- 
lated angular dependence of the two AFMR modes does 
not match the experiment [Fig.lHJd)]. The predicted non- 
monotonic dependence of the S2 mode is a consequence of 
a field- induced phase transition, occurring around 5.5 T 
[3^ for 50° < 9 < 90° [see the discontinuity of curves 
in Fig. IHc)[. The apparent absence of the transition 
at exactly 9 — 90° is accidental, which can be deduced 
from Fig. SJc). The transition corresponds to the field 
where the lowest branch reaches zero frequency. Such 
behavior was indeed observed and justified with the SI 
model on another Fe-based triangular lattice [3J|. The 
absence of any irregularities in low-temperature magne- 
tization isotherms speaks against such transition in BN- 
FSO at least up to 23 T Q. We stress that, on the 
other hand, no such excitation softening exists for the 
DM model [Fig. |4lja)[. Finally, the SI model predicts the 
selection of the {en, ea) = 1) GS independent of the 
structural chirality. This contradicts with the experimen- 
tal observation of the (1, -1) GS in a crystal with struc- 
tural left handedness [18|, clearly making the SI model 
incompetent of selecting the proper ground state. This 
selection is independent of the structural chirality, be- 
cause it is a single-ion property. 

The DM model [Eq. ([S])[, on the other hand, yields ex- 
cellent agreement with the experiment for \dc\ =45 mK. 
Both the V — B diagram [Fig. [4Ua)[ and the finer-detail 
angular dependence of both AFMR modes [Fig.lH^d)] are 
simulated convincingly. The DM anisotropy explains 
the ESR line width and the finite gap to the Si branch. 
Its relative size \dc/ Ji \ — 0.4% is in good agreement with 
the order-of-magnitude prediction i^g/g = 0.25%, as de- 
rived by Moriya [1[. We note that the zero-field gap 
could also be explained by the presence of symmetric 
anisotropic exchange (AE) of a similar size \J^'^ / Ji\ — 
0.4 — 0.8% [2^. However, as the AE anisotropy re- 
sults from the second order perturbation in A (while 
DM results from the first order perturbation), it is of 
the size (1| {Ag/g)'^Ji = 6 • lO^^Ji and therefore neg- 
ligible. We stress that the in-plane DM component dab 
does not affect the detected lowest branches and can- 
not be determined from the studies of these excitations. 
This makes the above-presented X-band ESR investiga- 
tion, yielding dab/dc = 2.6, invaluable in determining 
also the \dab\ = 120 mK DM component. There is a very 
good agreement between the excitations observed in ESR 
and calculated with the DM model on one hand and the 



modes identified from the INS experiments [22|, |36| on 
the other hand. The calculated Sq branch [Fig. lU^a)] is 
the gapless Goldstone mode emerging at the origin of the 
reciprocal space (scattering vector Q=0) for the correla- 
tions of the c-axis spin components. The Si mode corre- 
sponds to the other, gapped excitations for chiral corre- 
lations of the a6-plane spin components at Q=0. These 
excitations at Q = ±q involve correlations of the c-axis 
spin-component [22| and cross Q=0 at « 400 GHz, yield- 
ing the S2 signal. 

The DM anisotropy thus proves to be the origin of 
the unique magnetic properties of BNFSO. Our simula- 
tions [25| further disclose that one of the two possible GS 
[Fig. Wlh)] is selected depending on the sign of dc- The 
experimentally detected GS (1,-1) is realized for dc > 0, 
when the DM energy [Eq. ^\ is minimized. Formation 
of a uniform vector chiral state at finite temperature is, 
however, challenged by formation of domain walls, which 
prosper on frustrated lattices due to thermal fiuctuations 
|37| . It is likely that these are responsible for diffuse 
neutron scattering coexisting with LRO in BNFSO far 
below T/v [3^. It is therefore striking that the minute 
dc — 45 mK term effectively imposes a single-domain GS 
at two orders of magnitude larger temperatures, as ex- 
perimentally verified at 1.5 K [l8|, and a chiral unbalance 
up to Tat, with chiral correlations persisting even in the 
paramagnetic state [2^ . 

In conclusion, we have shown that the selection of the 
unique chiral magnetic ground state in BNFSO is due to 
the Dzyaloshinsky-Moriya interaction as the dominant 
source of magnetic anisotropy. Although this term is 
small {dc = 45 mK, dab = 120 mK), it effectively over- 
comes thermal fluctuations and leads to a fully chiral 
state in structurally enantiopure crystals, whereas spin- 
liquid signatures might also be an ticip ated above T/v, 
possibly with skyrmion texturation [39l, |40| . 

We thank S. de Brion for fruitful discussion and ac- 
knowledge the financial support of the Slovenian Re- 
search Agency (projects Jl-2118 and BI-US/09-12-040). 
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ESR ABOVE Tjv 



In the linear response theory an electron spin resonance 
(ESR) spectrum is given by the Fourier transform of the 
relaxation function [l| 



fit) 



{M+{t)M-{0)) 
{M+{0)M-{0))'' 



(4) 



where — gi^s'^j S^ are transverse magnetization 
operators and (...) denotes canonical averaging. Line 
position (g-value) and line width of the spectrum con- 
tain important information on local anisotropies . The 
principal axes of the 51-tensor are set by a local symmetry 
of the crystal field. In BNFSO one of the principal axes 
(labeled Z) is parallel to a two-fold rotational axis pass- 
ing through the center of each Fe04 tetrahedron, while 
the other two (labeled Y and X) are parallel with the 
two normals to the Oj-Fe-Oj planes, j — 2,3 labeling two 
pairs of oxygens at a distance of 1.909 A and 1.856 A from 
Fe, respectively. For site 2 (see Fig. [5] and Fig. 1 of the 
main text) the Z axis corresponds to the crystallographic 
a direction, while X and Y axes are tilted by 30° and 
120° from the b* axis in the crystallographic b*c plane. 
The three magnetically non-equivalent Fe sites on each 
triangle yield a single ESR line due to exchange narrow- 
ing. Therefore, an axially symmetric effective (7-tensor is 
obtained after averaging over the three sites. We find 



gab — 



gzzcos a + gyysm a, 

gxx + gYYCos^a + gzzsin^a 



(5) 
(6) 



and g{0) = glcos^O + g^^siT?0, where 6 denotes the 
polar angle between the crystallographic c axis and the 
applied magnetic field. 

When the magnetic anisotropy term T-L' of the spin 
Hamiltonian [see Eq. (1) in the main text] is inferior to 
the isotropic exchange term T-Le and/or the Zeeman term 
Hz it can be treated as a perturbation. In cases of mod- 
erate line widths (A_B ^ ksJ/gp^B) the well-established 



Kubo-Tomita approach applies for T > J. The ESR ab- 
sorption spectrum then has a Lorentzian line shape and 
a full- width-half-maximum (FWHM) line width, given at 
infinite temperatures as 
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The second and the fourth moment of the resonance spec- 
trum are given by [3| 



M2 - 
M4 



{[n',s+] [s-,n']) 

{S+S-) 

{[H - Hz, [W, S+]] [H - Hz, [W, S-]] 
{S+S-) 



(8) 
, (9) 



where [ , ] denotes a commutator. The moments given by 
the spin Hamiltonian are always finite, while those of the 
Lorentzian line shape are infinite. Therefore, the true ex- 
perimental line shape deviates from the Lorentzian shape 
and a modified line shape has to be employed. Typ- 
ically, this deviation becomes notable only in extreme 



(a) • 




FIG. 5: (color online). Triangle of magnetically non- 
equivalent Fe04 tetrahedra in BNFSO in (a) the ah and (b) 
the b*c crystallographic plane. The two crystallographically 
non-equivalent oxygen sites are labeled as 02 and 03. The 
principal axes of the g-factor on the Fe site 2 are denoted by 
X, Y and Z; the latter is parallel to a local two-fold rotational 
axis. 
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wings. The criterion for the selection of a suitable line 
shape should be the agreement between the theoretical 
exchange field — kg /guB \/ M4/M2 , determined by 
the spin Hamiltonian, and the one given by a particular 
line shape I^Si]. We find that in BNFSO the approximation 
where the Lorentzian line shape is multiplied with the 
broad Gaussian, G{B) (x exp[ — {B — Bq)^ /2i?e] , yields 
an excellent agreement of the exchange fields. This situ- 
ation is usual in systems with strong exchange, where the 
Gaussian time decay of spin correlations can be assumed. 



This approximation yields the •\/27r prefactor in Eq ([7]) . 

Although tedious, calculations of the two moments 
lead to analytic expressions of the ESR line width. We 
calculated it separately for the case of the single-ion (SI) 
anisotropy, taking into account spin operators up to the 
fourth order [see Eq. (2) in the main text], and for the 
Dzyaloshinsky-Moriya (DM) interaction [see Eq. (3) in 
the main text]. For the spin lattice shown in Fig. 1 of 
the main text we derived the following expressions for 
the ESR line width 



AB^ = V2^ 



AS, 



DM 



27 [AD-^ + fa^ 



10 p2 
9 



faF- 



2 (18D2 + 5F^ + 15aF) cos 20] ^ 



gfiB \ 20 J|j [189L»2 + 720a^ + 175F^ + 525aF - | (27i:>2 + 25F^ + 75aF) cos 29] ' 
I 105 [5d2^ + 6rf2 + (rf2^_2d2)cos 2ef 



gfiB V 32 [35 J2^d2, + ejg^dl + {2JS^dl - TJ^M cos 29] 



(10) 



fll) 



Here J|j = + 2 J| + J| + J| + J|, = 3 + 2 J| + 
J| + Jl + J| and Jg^ = 18 + 14 J| + 7J| + 7 J| + 7 J|. 

The ratios of the anisotropy parameters within each 
model are set by the temperature independent experi- 
mental line-width ratio A.B±/ A.B\\ — 1.27. In case of 
the SI model, our calculations predict AB^/ /S.B\\ = 1 for 
the a term and Ai3j^/Ai?|| ~ 1.5 for the D term. Due to 
its cubic symmetry the former term does not yield any 
anisotropy, however in still gives a finite ESR line width, 
because it does not commute with Hz- In case of the DM 
model, the dab anisotropy leads to AB±/ AB^^ — 1.5 and 
the dc anisotropy to AB±/ AB^^ — 0.5. In both models, 
a single dominant parameter thus fails to account for the 
experiment. 



AFMR BELOW Tat 

A mean-field approach is used to calculate the antifer- 
romagnetic resonance (AFMR) modes. In this approx- 
imation, sublattices are introduced instead of individu- 
ally treatment of the spins. The corresponding magne- 
tizations are given by Mj = —NgfiBi^j), where N is 
the number of Fe'^''" magnetic ions in the j-th sublat- 
tice. The AFMR modes are calculated within the spin 
Hamiltonian given by Eqs. (l)-(3) in the main text. Ap- 
plying the mean-field approximation (3| we can rewrite 
this Hamiltonian into a magnetic free energy J- per Fe 
site. 



E 



m: 



i>] 



^ {35Mlz - Mlz{30M' - 25}) 



const. 



(12) 



where the sums run over all sublattices. The molecular 
field constants are defined as 



^ ij J ij 



N{giJ.B) 



2 ' 



N{g^iB 



— F- 



N{gfiB 



,4 ' "'J 



N{g^iB) 



(13) 



where Zij is the number of the ij neighbors. 



The magnetic ground state is obtained from minimiza- 
tion of the free energy. This then yields effective mag- 
netic fields acting on individual sublattice magnetization 
Bj = — , which is introduced into the equations of 
motion for all sublattices 



dt 



9^J■B 



(14) 
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FIG. 6: (color online). Dependence of the energy gap to 
the lowest gapped branch on the anisotropy size A for the 
Dzyaloshinsky-Moriya (DM), single- ion (SI) anisotropy and 
symmetric anisotropic exchange (AE). The horizontal line 
corresponds to the experimentally observed gap. 

To make the problem tractable, harmonic oscillations are 
assumed, dMj/dt cx exp (icut)^ and the equations are lin- 
earized. The problem is solved numerically - the param- 
eters Jij, dij, Dj, dj and Fj are adjusted to fit the field- 
and the angular dependence of the antiferromagnetic res- 
onance modes. 

In this mean-field approach all information about the 
interactions between spins within a particular sublattice 
is lost. Therefore, a meaningful choice of sublattices is of 
great importance. The ferro-triangular magnetic order 
allows a restriction to six sublattices in each triangular 
plane, in order to account for the lowest lying magnon 
modes (see Fig. 1 in the main text for the site number- 
ing) . Thus both the intratriangle exchange Ji and the in- 
tertriangle exchange J2 are taken into account. The effect 
of the latter is however simplified, since all the l'-3' tri- 
angles neighboring the basic 1-3 triangle are assumed the 
same, i.e., effectively, only one 1-3 and one l'-3' triangle 
with the inter-triangle exchange J2 are needed. Finally, 
to take into the account also the interlayer interactions 
( J3, J4 and J5), we stacked 14 pairs of 1-3 and l'-3' trian- 
gles and used periodic boundary conditions, thus allowing 
for the experimental pitch of 2tt/7 along the c-axis. 

The resulting magnetic ground state is a helix with 
the ordering vector q = (0, 0, 1/7). The plane of the he- 
lix will, however, be randomly oriented if H' is set to 
zero, thus resulting in a sinusoidal tilt (along c) of spins 
on a given site from the c axis. This remains the case if 
a finite in-plane DM component dab is present. Adding 
a finite out-of-plane DM component dc aligns the helix 
plane with the crystallographic ab plane. The triangular 
chirality is found anti-parallel to the dc vector, thus mini- 
mizing the DM interaction energy. The helicity is always 



anti-parallel to the triangular chirality due to the left- 
handedness of the crystal structure imposing J5 > J3. In 
the SI model, on the other hand, a > leaves the helix 
plane randomly oriented, while D > randomly imposes 
helix planes with their normals perpendicular to the c 
axis. On the other hand, spins are bound to ab planes 
for D < and/or a < 0, which, however, always select 
the state with positive triangular chirality (and negative 
helicity) . 

The AFMR modes that suit the experimental obser- 
vations best are plotted in Fig. 4 of the main text for 
both anisotropy models. In zero field, the three low- 
est modes define the lowest two branches Sq and Si, 
which are characterized by J2a ~ ^ ^^'^ ^ linear mo- 
tion of spins. The crucial difference between the two 
anisotropy models is that in the DM model the So mode 
is the gapless Goldstone mode and the two Si modes oc- 
cur at (A = 73 GHz). In the SI model, on the other 
hand, all three modes are gapped - a two-fold degener- 
ate mode is found at A = 61 GHz and the third mode at 
A = 66 GHz. The small splitting is due to the fact that 
the SI anisotropy in BNFSO is multi-axial (see Fig. 1 of 
the main text). If the magnetic helix was incommensu- 
rate a gapless Goldstone mode would exist also in the 
SI model, as there would be no energy cost in rotating 
the helix around its axis. The gapless mode is, however, 
absent for a commensurate helix. Experiments imply the 
latter. Namely, no changes in the propagation vector, 
which is a rational fraction of c with small numerator 
and denominator co-primes, were detected in the tem- 
perature range of the ordered phase 0|. This is not in 
line with the incommensurability, arising from a delicate 
balance between the exchange parameters J3 — J5, then 
subject to strong sensitivity on temperature @. This 
then justifies the use of the finite lattice of 14 sites along 
the c axis in our modeling. 

The zero-field gap to the lowest gapped mode opens in 
a square-root fashion (Fig. ^ for both the DM and the 
SI model. We have calculated the dependence of the gap 
also for symmetric anisotropic exchange (AE) 

n'AE = J2 ^E^f 5|, (15) 

(y) 

where the sum runs over the nearest-neighbor Fe pairs, 
and f either corresponds to crystallographic c direction 
or a general local direction with respect to the ij bond 
within the ab plane. Similarly to other two models, in 
the AE model the square-root dependence of the gap is 
observed. Moreover, all models predict a very similar size 
of the anisotropy A, which is required to explain the gap. 
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